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Abstract 

In recent work, Zeilberger and the author used a functional equations approach for enu- 
| merating permutations with r occurrences of the pattern 12 . . . k. In particular, the approach 

■ yielded a polynomial-time enumeration algorithm for any fixed r > 0. We extend that ap- 

proach to patterns of the form 12 . . . (ft— 2)(fc)(fc— 1) by deriving analogous functional equations 
and using them to develop similar algorithms that enumerate permutations with r occurrences 
of the pattern. We also generalize those techniques to handle patterns of the form 23 ... k 1 and 
derive analogous functional equations and enumeration algorithms. Finally, we show how the 
functional equations and algorithms can be modified to track inversions as well as handle mul- 
l~i ■ tiple patterns simultaneously. This paper is accompanied by Maple packages that implement 

' the algorithms described. 

1 Introduction 

Let a = (Ji . . . (Jfe be a sequence of k distinct positive integers. We define the reduction red(a) 
to be the length k permutation t = ti .. . Tfc that is order- isomorphic to a (i.e., cr, < <jj if and 
only if Tj < Tj for every i and j). Given a (permutation) pattern r G Sk, we say that a permu- 
, tation 7T = 7Ti . . . ir n contains the pattern r if there exists 1 < ii < 12 < . . . < ik < n such that 

red(7Ti 1 7Ti 2 . . . m k ) = r, in which case we call n^ir^ . . . ni k an occurrence of r. We will define N t (tt) 
to be the number of occurrences of r in 7T. For example, if the pattern r = 123, the permutation 
53412 avoids the pattern r (so iVi 2 3 (534 1 2) = 0), whereas the permutation 52134 contains two 
occurrences of r (so A^i23(52 1 34) = 2). 

For a pattern t and non-negative integer r > 0, we define the set 

<S n (r, r) := {ir G S n : ir has exactly r occurrences of the pattern t} 
and also define s n (r, r) := \S n (r,r)\. The corresponding generating function is defined as 

oo 

F r T {x) :=^ S „(T,r)/. 

n=0 

Observe that the classical pattern avoidance problem corresponds to the case where r = and has 
been well studied. In this setting, S n (r, 0) is known to be enumerated by the Catalan numbers if 
t E S3 [11U19) . However, much is still unknown even for length 4 patterns. For example, permuta- 
tions avoiding the pattern 1324 have been notoriously difficult to enumerate. Precise asymptotics 
are not even known, although Bona recently gave an improved upper bound for the growth rate 
(in [3]) by modifying the approach used by Claesson, Jehnek, and Steingrfmsson [8]. The survey 
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paper by Kitaev and Mansour |10j provides an extensive overview of work in this area as well as 
related problems in permutation patterns. 

While the more general problem (where r > 0) has also been studied, the work has usually been 
restricted to small patterns (usually length three) and small r. In [Tl], Noonan studied permu- 
tations containing exactly one occurrence of 123 and proved that s n (123,l) = f („ 2 ™ 3 )- Burstein 
recently gave a short combinatorial proof for the result [6]. In [15], Noonan and Zeilberger pre- 
sented an approach using functional equations to enumerate s n (r, r) for small r and for the patterns 
123, 312, and 1234. Subsequent work has been done by Bona 011], Fulmek [9], Mansour and Vain- 
shtein [12], Callan [7], and many others. Many of these focused on finding F£(x) for r £ S3 and 
for small r. 

One difficulty arising from the initial Noonan-Zcilberger functional equation approach in |15] 
was that the approach became very complicated for even r = 2. In addition, there are many 
patterns that this approach does not readily extend to. One such pattern (explicitly mentioned 
in [IS]) is 1432. A modified approach was recently presented in [13] for the case of increasing 
patterns. Given a fixed r > 0, the resulting enumeration algorithm for computing s„(12 ...k,r) 
is polynomial-time (in n). This, in a sense, tackles the first difficulty from [15] and allows us to 
enumerate the sequence s„(12 . . . k, r) for even larger fixed r. 

In this paper, we extend the enumeration techniques in |13) to new families of patterns (includ- 
ing the pattern 1432) as well as multiple patterns. It should be noted that this general approach 
is different from the enumeration schemes approach pioneered by Zeilberger |21] and extended by 
Vatter [20] , Pudwell [TH [2] , and Baxter [TJ [2] . The enumeration schemes approach is useful for 
enumerating pattern-avoiding permutations (the r = case) but does not appear to be readily 
adaptable to the generalized setting for permutations with r > occurrences of a pattern. 

The paper is organized in the following manner. Section [2] extends the approach in [13] to the 
patterns 132, 1243, 12354, and so on. Section [3] generalizes the techniques used in [13] and applies 
them to the patterns 231, 2341 (which is equivalent to 1432), and so on. Section [4] extends this 
approach to handle multiple patterns simultaneously as well as refining by the inversion numbed- 
Section [5] lists some possible future work as well as some conjectures. The enumeration algorithms 
developed in this paper are implemented in the Maple packages FINCR, FINCRT, F231, F2341, 
F123nl32, F1234nl243, and FS3. They are all available from the author's website. 

2 Counting occurrences of the pattern 12 ... (h — 2){k){k — 1) 

In this section, we adapt the approach in [13] (for increasing patterns) to the patterns 12 . . . (k — 
2)(k)(k — 1). We first handle the case of 132 in full detail and then outline how to generalize this 
approach to patterns 1243, 12354, and so on. 

2.1 Permutations containing 132 

Given a (fixed) pattern r and non-negative integer n, we define the polynomial 

/„(<):= J2 tNrM - (!) 

lr This is technically the same as tracking the number of 21 patterns that occur. 
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Observe that the coefficient of t r in f n {t) is exactly equal to s n (r, r). For a fixed pattern r and 
fixed r > 0, our goal is to quickly compute s n (r, r). In the remainder of this section, we will 
assume that t = 132. 



In addition to the variable t, we introduce the catalytic variables xx,...,x n and define the 
weight of a length n permutation ir = tti . . . n n to be 



weight 132 (7r) : =^W jj 



#{(et,i>) : 7r a >ir i) =i, l<o<6<n} 



In general, this will be written more simply as weight^) when the fixed pattern is clear from con- 
text (in this case 132). For example, weight(12345) = 1, weight(13245) = tx 2 , and weight(25143) = 

t X -yX r^X 

For each n, we define the polynomial 

P n (t;xi,...,x n ) := 2J weight(7r). 

Observe that P n is essentially a generalized multi-variate polynomial for /„ and in particular, 
P n {t;l,...,l) = f n (t). 

Given a permutation n = 7Ti . . . 7r„, suppose that tti = i. Note that N^i^) is equal to the 
number of occurrences of 132 in -k 2 . . . 7r„ plus the number of occurrences of 21 in -ki . . . 7r„, where 
the term corresponding to the "1" is larger than i. Letting n' be red(7T2 . . .7r„), we observe that 

weight(7r) = x%X2 ■ ■ • weight(7r')U < ->-tx i+ i,xi + i->-t : i: <+ 2,...,a„_i-j-tx„- (2) 
This leads to the functional equation: 

n 

P n {t\xi,. . . ,x n ) = ^xix 2 ■ ..Xi^i ■ P n -i(t;xi, . . .,Xi-i,tx i+1 ,. . .,tx n ). (FE132) 

Once P n (t\x\, . . . ,x n ) is computed, the catalytic variables x\, ...,x n can all be set to 1 to get 
f n (t) = P„(t;l,...,l). 

However, it is not necessary to compute P n (t; x\, . . . , x n ) in its entirety prior to setting the 
catalytic variables to 1. Observe that by (IFE132p . we have: 

n 

P n (t; 1, . . . , 1) = P n _i(£; 1 [i — 1 times], t [n — i times]). 

i=l 

We get terms of the form P ao+ai (t; 1 [do times], t [a± times]) in the summation, which can again 
be plugged into (|FE132j) to get: 

Pa a +a x {t\ 1 [o-o times], t [<Zi times]) = P ao+Ol _i(l [i — 1 times] ,t [ao — i times] ,t 2 [a\ times]) 

+ t t ~ 1 P ao+ai _i(l [ap times], t [i — 1 times], t 2 [a\ — i times]) 



i=l 



Now, we must deal with terms of the form P ao+ai+a2 (t; 1 [ao times], t [a\ times], t 2 [a 2 times]). 
We can continue this recursive process of plugging new terms into (|FE132p to eventually compute 
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fn{t) = P n (t; 1 [n times]). This is much faster than the direct weighted counting of all n\ permu- 
tations, although it is still unfortunately an exponential-time (and memory) algorithm. 

This algorithm has been implemented in the procedure F132f ull(n,t) (in the Maple package 
FINCRT). For example, the Maple call F132full(8,t) ; computes f 8 (t) and outputs: 

t 31 + 7 t 30 + 20 t 2S + g 7 t 27 + 41 t 26 + 1Q g f 25 + lg2 f 24 + lg g f 23 + 333 f 22 + 397 ; 21 + M? ; 20 + ygQ t 

+ 1048 i 18 + 1152 t 17 + 1417 t 16 + 1576 t 15 + 1770 1 14 + 1853 i 13 + 2321 i 12 + 2088 t 11 + 2620 t 10 
+2401 1 9 + 2682 i 8 + 2489 t 7 + 2858 t 6 + 2225 t 5 + 2593 i 4 + 1918 t 3 + 1950 t 2 + 1287 < + 1430 

Suppose that for a small fixed r > 0, we wanted the first 20 terms of the sequence s n (132,r). 
By this functional equation approach, one would compute f n (t) and extract the coefficient of t r 
for each n up to 20. This approach would expend quite a bit of computational effort in gener- 
ating unnecessary information (namely, all the t k terms where k > r). This issue can mostly be 
circumvented, however, by a couple of observations. The first is the following lemma from |13j : 



Lemma 1. Let n = ao + a\ + . . . + a s (where ai > for each i ) and suppose s > r + 1. Then, the 
coefficients of t , i 1 , . . . ,t r in 



This lemma allows us to collapse all the higher powers of t into the t r+1 coefficient and allows us 
to consider objects of the form P n (t; 1 [ao times] , . . . , t r [a r times] , t r+1 [a r +i times]) regardless of 
how large n is. 

Let n :— oq + ai + . . . + a r+ \. Also, for any expression R and positive integer k, let R$k denote 
R [k times]. For example, t 3 $4 is shorthand for t 3 ,t 3 ,t 3 ,t 3 . Now for any polynomial p(t) in the 
variable t, let p^ r \t) denote the polynomial of degree (at most) r obtained by discarding all powers 
of t larger than r. Also, define the operator CHOP r by CHOP r [p(i)] := p {r) (t). 

An application of (|FE132p and CHOP r to Pn^ (t; 1 [a times] , . . . , t r [a r times] , t r+1 [a r+1 times]) 



P n (t; 1 [ao times] , . . . , t s 1 [a s _i times] , t s [a s times]) 
P n (t, 1 [ao times] , . . . , t r [a r times] , t r+1 [a r+ i + a r+ 2 + . . . + a s times]) 



all vanish. 



Proof. The more general function P n {t;x\, . . . , 



x n ) is a mufti- variate polynomial. 



□ 



becomes: 



Pi r) {t; l$o , • ■ • , t r $a r , t r+1 $a r+1 ) 



= CHOP r p n-i(*; - 1). *S(<«o - i),t 2 $a u t r $a r _i, f +1 $(a r + a r+1 )) 



+ ^ t i-1 Pi!2i(t; l$ao, - 1), t 2 $(ai - i), i 3 $a 2 , • ■ • , t r $a r _i, f +1 $(a r + a r+ i)) 



02 



+ t ai+2 ^ P { n r \{t; l$a , t$oi, t 2 ${i - 1), t 3 $(a 2 -»),.-., t r $a r -i, t r+1 $(a r + a r+ i)) 



+ 



+ J2 f »i+2a,+...+ra r +(r+l)(i-l)pM i ^. 1$aQ) ^ ... ; r +1 $(a r+1 - 1)) . 



i=l 
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Due to the CHOP,- operator, many terms automatically disappear because of the power of t 
in front. From a computational perspective, this observation eliminates many unnecessary terms 
and hence circumvents a lot of unnecessary computation. This has been automated in the Maple 
package FINCRT so that a computer can derive a "scheme" for any fixed r (completely on its own) 
and use it to enumerate s n (132,r) for as many terms as the user wants0 

For example, the Maple call F132rN(5,15) ; for the first 15 terms of s„(132, 5) produces the 
sequence: 

0, 0, 0, 0, 5, 55, 394, 2225, 11539, 57064, 273612, 1283621, 5924924, 27005978, 121861262 



2.2 Extending to the pattern 1243 

First, we outline how to extend the previous approach to the pattern 1243. In addition to the 
variable t, we now introduce 2n catalytic variables x\, . . . ,x n and y\, . . . , y n . The weight of a 
length n permutation tt = w± . . . n n will now be 

n 

Weight(7r) ■= t Nl2i3 ^ JJx*^"'^ 1 7ra>nb=i ' ^-<1<b<n} , y#{( a - b , c ) ■ n a =i<TT a <n h , l<a<b<c<n} 
i=l 

For example, weight(123456) = 1 and weight(135624) = t 2 x\x\y\y\. 

For each n, we define the polynomial 

Pn(t;x!, ...,x n ,yi,.- .,y n ) ■= ^ weight(Tr). 

7r£S„ 

Now given a permutation n = ~k\ . . . 7r n , suppose that w± — i. Observe that iVi243(7r) is equal to 
the number of occurrences of 1243 in 7T2 . . . 7r„ plus the number of occurrences of 132 in 7T2 . . . n n 
where the term corresponding to the "1" is larger than i. Letting ir' := red(7T2 . . . 7r„), we observe 
that 

weight(7r) = x 1 x 2 ...x i -i ■ weight(Tr')\ Xi ^ yiXi+lt ... !Xn _ 1 ^y iXn! y i ^ty i+1 ,...,y n - 1 ^tv n - (3) 
This leads to the functional equation: 

P n (i, X\ , . . . , x n , yi , . . . , yn) — 

^ xix 2 ■ ■ ■ Xi-i ■ P„-i(t; xx, . . . , Xi-x,yiX i+ i,. . . , y l x n , yi, . . . , y^-i, . . . , ty n ). (FE1243) 

i=l 

Again, our goal is to compute f n {t) = P n (t; 1 [2n times]). We can apply the same computational 
methods as before. For example, we can apply (|FE1243p directly to P ra (i; 1 [2n times]) (and more 
generally, to objects of the form P n (t; 1 [do times] , . . . , t Sl [a Sl times] , 1 [bo times] , . . . , t S2 [b S2 times])) 
to compute f n (t) . This again gives us an algorithm that is faster than the direct weighted counting 
of n! permutations but is still exponential-time (and memory). 

This algorithm has been implemented in the procedure F1243full(n,t) (in Maple package 
FINCRT). For example, the Maple call F1243full(8,t) ; computes f s (t) and outputs: 

f 36 + t 31 + 10 i 30 + 3 1 2& + 13 1 27 + 9 1 26 + 8 1 25 + 37 i 24 + 16 1 23 + 16 < 22 + 49 i 21 + 60 1 20 
+41 i 19 + 130 1 18 + 81 1 17 + 157 1 16 + 266 1 15 + 184 1 14 + 233 1 13 + 542 i 12 + 356 1 11 + 771 i 10 
+877 1 9 + 975 t s + 972 1 7 + 2180 1 6 + 1710 1 5 + 2658 1 4 + 3119 1 3 + 4600 1 2 + 4478 * + 15767 



2 The "scheme" mentioned here is a liberal application of the word and differs from enumeration schemes. 
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Additionally, both the obvious analog of Lemma[T]as we U as the computational reduction using 
the CHOP r operator still apply in this setting. This has also been automated in the Maple package 
FINCRT. 

For example, the Maple call F1243rN(l , 15) ; for the first 15 terms of s n (1243, 1) produces the 
sequence: 

0, 0, 0, 1, 11, 88, 638, 4478, 31199, 218033, 1535207, 10910759, 78310579, 567588264, 4152765025 
and the Maple call F1243rN(2, 15) ; for the first 15 terms of s„(1243, 2) produces the sequence: 
0, 0, 0, 0, 4, 56, 543, 4600, 36691, 284370, 2174352, 16533360, 125572259, 955035260, 7283925999 

2.3 Extending to longer patterns 

The approach for the patterns 132 and 1243 can be extended analogously to longer patterns of 
the form 12 ... (k — 2)(k)(k — 1). For example, if the pattern r = 12354, we consider the variable 
t and 3n catalytic variables: xi, . . . ,x n and t/i, . . . , y n and z±, . . . , z n . The weight of a length n 
permutation n = tti . . . ir n will now be 

n 

Weight(7r) = t Nl235i ^ J^X # ^ a:6 ' 1 * a>7Tb=i } . y#{( a > & ' c ) : Ta=i<Te<W!.} , z #{(a,b,c,d) : ir a =i<T! b <Tt d <Tt c } 
i=l 

where it is always assumed that a < b < c < d. 

An analogous functional equation is derived for the corresponding polynomial 

) := ^ weight (tt) 

lt£5„ 

and all the analogous computational methods work in this setting as well. The 12354 case has also 
been automated in the Maple package FINCRT. 

For example, the Maple call F12354rN(0, 14) ; for the first 14 terms of s„(1243,0) produces 
the sequence: 

1, 2, 6, 24, 119, 694, 4582, 33324, 261808, 2190688, 19318688, 178108704, 1705985883, 16891621166 
and the Maple call F12354rN(l , 15) ; for the first 15 terms of s„(1243, 1) produces the sequence: 
0, 0, 0, 0, 1, 19, 246, 2767, 29384, 305646, 3170684, 33104118, 349462727, 3738073247, 40549242195 

3 Counting occurrences of the pattern 23 . . . kl 

In this section, we extend the previous techniques to handle patterns of the form 23 . . .kl. Although 
s„(231,r) — s ra (132,r) for every r and n (by reversal), we will develop an approach for handling 
231 directhjfl and then show how this can be extended to longer patterns of the form 23 . . . fcl. This 
new approach for handling 231 will also be necessary in the next section for considering multiple 
patterns simultaneously. 



As opposed to computing the equivalent pattern 132. 
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3.1 Permutations containing 231 

In this section, we will assume that our (fixed) pattern r = 231. We define the analogous polyno- 
mial 



f n (t) := ]T t N ^ 



Recall that the coefficient of t r in f n {t) will be exactly s„(231,r). 

In addition to the variable t, we introduce n(n+ 1)/2 catalytic variables x t j with 1 < j < i < n 
and define the weight of a permutation tt = m . . . 7r„ to be 

#{(a,b) : 7r a >7Tf,, 7r a — i, 7Tb<j, l<a<6<n} 



n 



weight(Tr) ^i^'W 

l<j<i<n 

For example, weight(12345) = 1 and weight(24153) = i 2 i£2. 22:4. 2^4, 32:4 ^5^X5 



We will again define an analogous multi-variate polynomial P n on all the previously defined 
variables. However, for notational convenience, the Xij variables will be written as a matrix of 
variables: 



X n :— 



X\A 



xi, 



X'i/i, 



XnA 



(4) 



where we will disregard the entries above the diagonal (i.e., the Xij entries where j > i). 
For each n, we now define the polynomial 

P n {t-X n ) := weight(Tr). 

Recall that P n (t; 1) = f n (t), where 1 is the matrix of all l's. 

Next, we derive an equation (analogous to Eq. [2] and [3} for this weight function. Given a 
permutation n = 7Tx . . . 7r n , suppose that m — i. Observe that N23i(ir) is equal to the number of 
occurrences of 231 in 7T2 . . . 7r n plus the number of occurrences of 21 in it2 ■ ■ ■ 7r„, where the term 
corresponding to the "2" is greater than i and the term corresponding to the "1" is less than i. 
We can also make the following two observations. First, in weight^), the exponents of Xk,i and 
Xk,i+i are equal for each k (since 7Ti = i). Second, the number of 231 patterns that include the 
first term tti = i is the exponent of x i+ i ti plus the exponent of Xi + 2,i plus . . . plus the exponent 
x n< i. Now, let 7r' be red^ . . . 7r n ), and define the collection of substitutions A to be 



Xb, c Xb+i,c b>i,c<i 
A := < x b ^ c -4 a;;, + i !C+ i b>i,c>i 
Xb, c -> txb+i, c ■ Xb+i,c+i b>i,c = i 



(5) 



and observe that 



weight (?r) 



r° r 1 



weight (77') j 



A- 



(6) 
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Now, define the operator R± on a square matrix X n and i < n to be: 



£1-1,1 
z*+i,i 



^i— 1 

£«+l,j-l ^i+l,i+2 



+ 1 



#n,i+2 



^l,n 



^-i — l,n 
<£i+l,n 



In essence, the i?i operator deletes the i-th row, merges the i-th and (i + l)-th columns via term- 
by-term multiplication, and multiplies this new column by a factor of t. If i = n, then Ri(X n ,i) is 
defined to be the (n— 1) x (n— 1) matrix obtained by deleting the n-th row and n-th column from X n . 



Observe now that the weight equation (Eq. [5]) leads to the functional equation: 

n 

P n (t; X n )=J2 4,1x1,2 •• -xtf -P n -i(t;Ri(X n ,i)). (FE231) 
i=i 

Note that while all entries in the matrix are changed for consistency, we will continue to disregard 
the entries above the diagonal. 

Again, our goal is to compute P n (t;l), and the analogous computational techniques from 
previous sections will also apply in this setting. For example, we can apply (|FE231[) directly to 
P n (t; 1) as opposed to computing P n (t; X n ) symbolically and substituting Xij — 1 at the end. The 
following result, which is obvious from the definition of the operator R\, provides a substantial 
simplification: 

Lemma 2. Let A be a square matrix where every row is identical (i.e., the i-th row and the j-th 
row are equal for every Then, R\(A,i) will also be a square matrix with identical rows. 

By Lemma [2] repeated applications of R\ to the all ones matrix 1 will still result in a matrix 
with identical rows. Therefore, it is sufficient to keep track of only one row as opposed to the 
entire matrix. Also observe that repeated applications of Ri to the matrix 1 will always result in a 
matrix whose entries are powers of t. Let Q n (t; c\, . . . , c n ) denote the polynomial P n (t; C), where 
C is the n x n matrix where every row is [c\, . . . , c n ] and every Cj is a power of t. This leads to a 
functional equation analogous to (|FE231[) : 

Q n (t;ci, . . . ,c„) = ^c?C2 . . .c l r x ■ Q n _ 1 (t;a, . . . , Ci_i, tCfCj+i, c i+2 , ■ ■ ■ ,c„). (FE231c) 

i=i 

Note that Q n (t; 1 [n times]) is exactly our desired polynomial P n (t; 1) = f n (t). However, this in- 
terpretation only forces us to deal with n catalytic variables (the Cj's) as opposed to n(n + l)/2 
catalytic variables (the a^j's). Just as in prior sections, we can repeatedly apply our functional 
equation (|FE231cl) to compute Q n (t; 1 [n times]). 

When the sequence s„(231, r) is desired for a fixed r, the obvious analog of Lemma[T]and the 
computational reduction using the CHOP r operator can again be used. This has been implemented 
in the Maple package F2310 

4 Although all output would be equivalent to the 132 case, the approach here will be necessary when considering 
multiple patterns. 
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3.2 Extending to the pattern 2341 

In this section, we outline how to extend the approach for 231 to an analogous (but more com- 
plicated) approach for 2341. In addition to the variable t, we now introduce n(n + l)/2 catalytic 
variables x^j with 1 < j < i < n and n(n + l)/2 more catalytic variables y^j with 1 < j ' < i < n 
(a total of n(n + 1) catalytic variables). Define the weight of a permutation tt — tt\ . . . ir n to be 



t 



N 2 3il (t) 



n 



weight (tt) :— 

#{(a.b) : 7r a >7Tb, 7r a — i, 7Tb<j, l<a<6<n} #{(a,b,c) : 7r c <7r a <7Tb, 7r a —i, 7r c <j, l<a<b<c<n} 



l<j<z<n 



For example, weight(24351) = i 2 ^, 2^3, 2^3, 3^4, 2X4,3a;i,4^5 : 2a;5,3a;5, 4^5, 52/1,22/3, 22/3,32/4,22/4, 32/4,4- 
The sc»j variables and the j/jj variables will be written as matrices of variables: 

2/i,i ••• 2/l,n 



5Cl,l 



X-n, 1 



2M,i 



(8) 



where we will disregard the entries above the diagonal. 



For each n, we define the polynomial 

P n (t;X n ,Y n ) := ^ weight(Tr) 

7r£S„ 

and again P„(i; 1, 1) = f n (t) is our desired polynomial. 

Next, we derive an equation (analogous to Eq. [2j [3j and |6]) for this weight function. Given a 
permutation tt = 7ri...7r n , suppose that 7Ti = i. Observe that -/V2341 (it) is equal to the number 
of occurrences of 2341 in TT2 ■ ■ ■ 7r„ plus the number of occurrences of 231 in TT2 ■ ■ ■ TT n , where the 
term corresponding to the "2" is greater than i and the term corresponding to the "1" is less than i. 

We can also make the following observations. First, in weight(7r), the exponents of Xk,i and 
Xk,i+i are equal and the exponents of yt,i and 2//M+1 are equal for each k (since tt\ — i). Second, 
the number of 2341 patterns that include the first term tt\ — i is the exponent of 2/i+i.i phis the 
exponent of 2/1+2, i phis . . . plus the exponent y n j. Third, the number of 231 patterns that include 
the first term tt\ = i (i.e., the "2" is equal to i) and whose "1" term is less than k is equal to the 
exponent of Xi+i.fc plus the exponent of Xi+2,k plus . . . plus the exponent of x n> }-. 

Now, let tt' be red(7T2 . . . 7r„), and define the collection of substitutions A' to be 



A' := I 



Xb, c -» yi,c • £&+i, c b > 

Xb,c -> 2)6+1. c+1 b > 
x b,c — > Vi,i ■ ^6+l,c • Xb+l.c+l 

Vb,c -> Vb+i.c b > 

2/b,c -> Vb+i.c+i b > 

V Vb.c -> ty b+ i.c ■ 2/b+i,c+i b>i,c 



1, c < 1 
i,c > i 
b > i, c — i 
,c<i 
,c> i 



(9) 
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and observe that 

weight(Tr) = xl lX \ 2 . . . x]- 1 ■ weight(7r')U<. (10) 

In addition to the previous R\ operator defined in Eq. [71 we define another operator R2 on two 
square matrices X n and Y n (of equal dimension) and i < n to be: 

X\ 1 ■ ■ ■ Xxi-X yi,iXX,iXx,i+X Xx,i+2 • ■ • Xx,n 



Xi—X,X Xi— X,i— X ' ' ' X$—X, 

R2(X n , Y n , i) :— Vi.l x i+XA '•" Vi,i— XXi+X,i-l Vi,iXi+X,iXi+X,i+X ^i+l,i+2 '•• ^i+l,n 



Ui : lXn,X ' ' ' yi,i — XXn,i — X ]Ji,iXn,iXn,i+X *^n,i+2 ' ' ' *^n,n 

(11) 

If i = n, then R,2(X n , Y n ,i) is defined to be the (n — 1) x (n — 1) matrix obtained by deleting the 
n-th row and n-th column from X n . 

Observe that this leads to the functional equation: 

n 

P n (t; X n , Y n ) = <i^,2 ■ ■ • x\f ■ P n -x(t; R 2 (X n , Y n ,i), Rx(Y n , »)). (FE2341) 

i=X 

As in prior sections, we recursively apply the functional equation directly to P n (t; 1, 1) (and sub- 
sequent instances of Pk). Observe that in this scenario, Lemma [2] still applies for the i?i operator 
and more specifically the "Y n " matrix in P n . While the lemma does not apply to the R2 operator, 
this still allows us to reduce the number of catalytic variables. Let Q n (t; C,dx, ■ • ■ ,d n ) denote the 
polynomial P n (t; C, D) where every entry of the nx n matrices C and D are powers of t and every 
row in D is [dx , ■ ■ ■ , d n ] . We derive an analogous functional equation: 

n 

Q n (t;C,dx, ...,d n ) = J~] c° ]C* i2 . . .c^ 1 ■ Q n -x(t; Rz{C, D, i), d x , . . . , c^-i, tdidk+x, ck +2 , . . .,d n ). 



(FE2341c) 



Using this recurrence to compute Q n (t; 1, 1 [n times]) will yield the desired polynomial f n {t). This 
approach allows us to deal with n(n + l)/2 + n catalytic variables (as opposed to n(n + 1) such 
variables). 

Additionally, for a fixed r, the sequence s„(2341, r) can be computed by applying LemmaUJand 
the CHOP r operator as necessary. This has been implemented in the procedure F2341rN(r,N) (in 
the Maple package F2341). 

For example, the Maple call F2341rN(l , 15) ; for the first 15 terms of s„(2341, 1) produces the 
sequence: 

0, 0, 0, 1, 11, 87, 625, 4378, 30671, 216883, 1552588, 11257405, 82635707, 613600423, 4604595573 
and the Maple call F2341rN(2, 15) ; for the first 15 terms of s n (2341,2) produces the sequence: 
0, 0, 0, 0, 5, 68, 626, 5038, 38541, 289785, 2172387, 16339840, 123650958, 942437531, 7236542705 
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While we do not present the details here, the same methodology can be applied to longer 
patterns of the form 23 . . . kl. Analogous functional equations can be derived and used for enu- 
meration. 

4 Further extensions 

4.1 Tracking inversions 

One of the most commonly studied permutation statistic is the inversion number. The inversion 
number of a permutation n = tti . . . 7r„, denoted by inv(7r), is the number of pairs such that 
1 < z < j < n and ~ki > itj. Equivalently, it is the number of occurrences of the pattern 21 in it. 
For a (fixed) pattern r, define the polynomial 

9n (t,q)~ J2 q invM t N ^l (12) 

Observe that g n (t, 1) is exactly f n (t) from before. 

Given a permutation ir = tt\ . . . n n , suppose that -k\ — i. Then, inv(7r) is equal to the number 
of inversions in 7Ta • . . ^ n phis the number of elements in 7T2, . . . , n n that are less than i. For any 
previously defined functional equation, it is enough to insert a factor in the summation. 

For example, if the fixed pattern is r = 132, the polynomial P n can be analogously defined as 

P n (t,q; Xl ,...,x n ) := ]T q invM ■ weight 132 (tt) 

and the analog to functional equation (|FE132[) would be 

n 

P n {t, q;xi,..., x n ) = 2J q l ~ x x\X2 . . . Xi-i ■ P n -i(t, <?; xi, . . . , Xi-i, tx i+ i, tx n ). 

i=l 

Similarly, the analogous functional equation to (|FE231c|) would be 

n 

Qn(t,q; ci, . . . , c„) = ^<7 i_1 c?C2 . . . c*" 1 • Q n _i(t, q; ci, . . . , Cj_i, tc t c i+1 , c t+2 , . . . ,c n ). 

i=l 

From here, all the previous computational techniques for quick enumeration still apply. 

This has been implemented in the procedures qF123r (n,r ,t ,q) and qF1234r(n,r,t,q) (in 
Maple package FINCR), qF132r (n,r ,t ,q) and qF1243r (n,r ,t ,q) (in Maple package FINCRT), 
qF231r(n,r,t,q) (in Maple package F231), and qF2341r (n,r ,t ,q) (in Maple package F2341). 

4.2 Counting multiple patterns in permutations 

In the preceding sections, various functional equations were derived by considering the first term 
of a typical permutation^! and deriving a recurrence. Given any collection of patterns where such 
recurrences can be derived, we can also consider those patterns simultaneously. 



5 As opposed to looking at the last term (or the smallest /largest term). 
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As an example, consider the patterns a — 123 and t = 132. The case of only the pattern 123 
was done in |13j . In this setting, the weight of a permutation tt = tt\ . . . 7r„ is defined by 



weight 123 (7r) := t NlwM JJ x\ 



._ fN 12 :,(Tv) TT ?? #{(<».6) : 7ra=i<7rf,, l<a<6<n} 
i=l 

and the corresponding polynomial is 



P n (t;x!, . . . ,x n ) := ^2 weight 123 (7r). 



ttG5„ 

The corresponding functional equation (referred to as the Noonan- Zeilberger Functional Equation) 
is 

n 

P n (t;xi, ...,x n )= y^a;? -i P n _ift;xi,.. .,Xi_i,ti i+ i, . . . ,te„). (NZFE) 

»=i 

This can be merged with the analogous quantities for 132 as follows. Let s and t be the variables 
corresponding to 123 and 132, respectively. Let x%, . . . , x n and y±, . . . , y n be two sets of catalytic 
variables, and define the weight of a permutation 7r = 7Ti . . . ir n by 

n 

Weight(7r) ■= s N ^3(T) i -N 13 2M a ,#l( tt . b ) : K a =i<-K b , l<a<b<n} y#{(a,b) : Tt a >-rt b =i, l<a<b<n} 

i=l 

For each n, we define the polynomial 

P n (s,t;xi, . . . ,x n ,yi, . . . ,y n ) := ^ weight(7r) 

7re5„ 

and can similarly derive the functional equation 

Pn j t : X \ , . . . , x n , yi , . . . , y^i) 

n 

^ x"~ l yxy 2 . . . j/i_iP„_i(s, t; n, . . . , a;<_i, 8X i+ i, . . . , sx n , yi,..., yi-i,ty i+1 , ty n ). 

i=l 

The same computational techniques from the prior sections apply here as well. 



This has been implemented in the Maple package F123nl32. For example, the Maple call 
F123rl32sN(2,2, 15) ; gives the first 15 terms of the sequence enumerating permutations with 2 
occurrences of 123 and 2 occurrences of 132: 

0, 0, 0, 1, 6, 26, 94, 306, 934, 2732, 7752, 21488, 58432, 156288, 411904 

and the Maple call F123rl32sN(4,2, 15) ; gives the first 15 terms of the sequence enumerating 
permutations with 4 occurrences of 123 and 2 occurrences of 132: 

0, 0, 0, 0, 1, 5, 23, 106, 450, 1740, 6214, 20831, 66427, 203550, 603920 

Other pairs (or larger sets) of patterns follow similarly, and the analogous 1234 and 1243 case 
has been implemented in the Maple package F1234nl243. 



Finally, it is possible to consider all length 3 patterns simultaneously. Only the patterns 123, 
132, and 231 were done directly, but analogous functional equations can be derived for 321, 312, 
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and 213. These six functional equations can be combined to count occurrences of all the length 
3 patterns. This has been implemented in the Maple package FS3. For example, the Maple call 
FS3full(7, [t[l] ,t[2] ,t[3] ,t[4] ,t[5] ,t[6]]) ; would produce the polynomial 



5 Conclusion 

In this work, we extended and generalized the techniques of [13] to the pattern families 12 . . . (k — 
2)(k)(k — 1) and 23 . . . fcl. In addition, we showed how this approach could be further extended to 
handle inversions and more generally, multiple patterns simultaneously. It would be interesting to 
see what additional patterns this approach can be applied to. Also, while the main results of this 
paper are enumeration algorithms, they are based off of rigorously derived functional equations. 
It would be interesting to find out if any additional information can be extracted from these func- 
tional equations. 

Finally, the techniques of this paper allow us to compute many sequences that lead to new con- 
jectures. Let c ryS {n) be the number of length n permutations with exactly r occurrences of 1234 
and exactly s occurrences of 1243. We will denote this as c(n) when r and s is clear from context. 
When r = s = 0, it is known that c(n) is exactly the Schroder numbers. For fixed r, s < 1, we are 
almost certain that c(n) is P-recursivqj (the algorithm can compute enough terms to guess a recur- 
rence). On the other hand, for the single pattern case, it is not clear if s„(1234, 1) is P-recursive |13j . 

Based off of empirical evidence, we also believe the following to be true: 

Conjecture 1. For each fixed r > and s > 0, the sequence enumerating length n permutations 
with exactly r occurrences of 123 and s occurrences of 132 is given by p(n)2 n , where p(n) is some 
polynomial of degree r + s. 

There are a number of results considering this type of problem (beginning with [17| 118]). but most 
such results limit themselves to s = 0, 1. If this general form were shown to hold for arbitrary r 
and s, the F123nl32 package could quickly compute enough terms to find explicit formulas and 
generating functions. 
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